Electromagnetic Shocks in Strong Magnetic Fields 
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We examine the propagation of electromagnetic radiation through a strong magnetic field using 
the method of characteristics. Owing to nonlinear effects associated with vacuum polarization, 
such waves can develop discontinuities analogous to hydrodynamical shocks. We derive shock jump 
conditions and discuss the physical nature of these non-linear waves 
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I. INTRODUCTION 

00 

■ The nonlinear properties of electromagnetic waves traveling through a magnetized vacuum is of particular interest 
in the study of neutron stars. Heisenberg and Euler (lj and Weisskopf j!| first derived nonlinear corrections to the 
Maxwell equations of the electromagnetic field. Lutzky and Toll H and Zheleznyakov and Fabrikant jf| applied the 
weak-field expansion to show that shock waves can develop in the electromagnetic field. Bialynicka-Birula || used 
the full expression for the nonlinear correction to the Lagrangian to study the generation of harmonics and other 
' nonlinear phenomena in the propagation of EM radiation. 

In this paper, we use the method of characteristics to study the evolution of waves governed by an arbitrary 
£NJ ■ Lagrangian and then apply these techniques to the Heisenberg-Euler-Weisskopf-Schwinger Lagrangian We 
^ \ find a concise expression for the opacity to shocking for disturbances travelling through an arbitrarily strong mangetic 
field. After deriving a form for the energy-momentum tensor for an arbitrary electromagnetic Lagrangian, we use 
relativitistic fluid mechanics to derive the shock jump conditions and the long-term evolution of electromagnetic 
. disturbances propagating in magnetic fields. 

o 

OO 

Q\ ■ II. DERIVING THE CHARACTERISTICS 

We will use the method of characteristics to study the evolution of a disturbance of the electromagnetic field. In 
general, the relativistic Lagrangian (£) of the electromagnetic field is a function of the two invariants of the field. We 



follow the notation of Lutzky and Toll H and Heisenberg and Euler and define 



/ = F^F^ = 2(|B| 2 -|E| 2 ) (1) 



and 
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. P. ; K = - \-e xp ^F Xp F^\ = -(4E • B) 2 . (2) 

As illustrated in Fig. |l|, we choose coordinates so that the radiation is polarized in the z-direction and travels along 
the y-axis toward positive y. The ambient magnetic field makes an angle <p with the electric field, and the projection 
of the magnetic field into x — y plane makes an angle 6 with respect to the ir-axis (magnetic field of the wave). 
With these definitions, the invariants are 

/ = 2 ((B cos 6 sin </> + B) 2 + (B sin 9 sin 0) 2 + (B cos </>) 2 - E 2 ) (3) 
K = -(4EB cos (f>) 2 (4) 

where B is the strength of the ambient magnetic field, E and B are the strengths of the electric and magnetic fields 
associated with the radiation. 

We characterize the traveling wave by a vector potential with one non-zero component, 
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A z = ip{y 7 t). 



(5) 



We have assumed that the vector potential has only one independent component which allows us to treat the problem 
using characteristics. Unfortunately, we cannot follow the interaction between two polarizations which was treated 
by Bialynicka-Birula JFJ. 

Using the coordinates, y and t, Hamilton's principle assumes the form 
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= 0. 
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where 



tp y = B and ipt 



-E, 



(7) 



and we have taken c = 1. 

Since the Lagrangian is a function of / and K alone, we can rewrite this in terms of derivatives of C with respect 
to / and K , / and K with respect to ip v and ipf and finally ip yy , ip y t, and i^tt by successive applications of the chain 
rule. Taking the partial derivatives yields an equation of the form, 
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This equation may be factored yielding 
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where 
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and a± is the speed of the travelling wave. 

We define two new coordinates (it, v) such that 



y = u + v and t = —{t + u + r_w) 



We now define the Ricmann invariants, 

V>+ = tp y - r+Vt 

1p-=tpy- T-1p t 

which are constant along curves with 
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To obtain a cursory understanding of the characteristics of this equation, we expand each of the coefficients about 
E, B = to first order yielding, 



a = a + aBipy - astpt + 0(B 2 ) 
b = b Q + b B ^ y -b E ^ t + 0(B 2 ) 
c = c + c B ip y - c E ipt + 0(B 2 ). 



The coefficients are 
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and 
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Using this linearization we estimate the magnitude of the inhomogeneous term in Eq. U and Eq. if we assume 
E t ~E and B t ~ B 



d__ d 

dy T± 



dt 



0(B 2 



so we can neglect it to first order. 
We estimate the Riemann invariants, 



V'+ = Aj - r+A + 0(B 2 ) 
= B + t+E + 0{B 2 ) 
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Now taking the limit where E and B themselves may be neglected we have 



o± = ±\/-a /c = ± 
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where 
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We can see from this result that r + = — t_ + O(B), therefore we have for the two Riemann invariants 

ip± = B ± r+£ + 0{B 2 ) 



(37) 



(38) 



Fig. H depicts how two adjacent characteristic may intersect. Additionally, we see that -0_ characteristics {i.e. lines 
along which is constant) that originate from regions without wave fields cross the ip+ characteristics. Therefore, 
we can argue that ip- = B — t + E + 0(B 2 ) = throughout the region to the right of the antenna because the fields 
are zero in this region. We can use the same argument for the region to the left of the antenna and find that in 
general B = t±E + 0(B 2 ) along the if>± characteristics. Furthermore, to first order, both ip± are constant along the 
characteristics; therefore, the slopes of the characteristics which depend only on ip+,ip- and the constant background 
field must be constants and the characteristics travel at a constant speed. 

Using the figure as a guide, we estimate the distance over which two adjacent characteristics can travel before 
intersecting is given by 
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Here we have used that fact that ip+ characteristics travel with velocity — r_ = t + + 0(B). 
To work with this equation further, we define an opacity due to shock formation and use B 
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where we have defined 



Bn = B cos 6 sin 6 and B r = B cos < 



(43) 



Where two of these characteristic curves intersect, a shock will form. To move further, we must choose the proper 
Lagrangian. 



III. THE NON-LINEAR LAGRANGIAN 



Heisenberg and Euler yj] and Weisskopf [pf independently derived the effective Lagrangian of the electromagnetic 
field using electron-hole theory. Schwinger [p[ later rederived the same result using quantum electrodynamics. In 
rationalized electromagnetic units, the Lagrangian is given by 
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In the weak field limit Heisenberg and Euler give 
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We define a dimensionless parameter £ to characterize the field strength (/) 

and use the analytic expansion of this Lagrangian for small K derived by Heyl and Hernquist : 
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The first two terms of this expansion are given by 
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where £W(aj) is the first derivative of the Riemann Zeta function. 

Using these definitions we can derive the various partial derivatives important for shock formation 



dC 
~dl 

dC 
dK 

d 2 C 
~dP 

d 2 C 



e 

he 



1 

~4 

^ P -E- 
288 he k 

1 e 2 



:X n 



16 he 



El 



(i) 



(1) 



r 1 



LS.V,, ( - 



6^ 



(2) 



9^o 



(i) 



fir- 



(2) 



1 e 2 



dldK 384 ftc 



i2.v (! ( - 



4X 



(2) 



- ) s r 4 - ( axf 



2X, 



(3) 



3X 



(2) 



d 3 £ 

a/ 3 



64 /ic fe 



3X 



IV. THE OPACITY TO SHOCKING 



Using the results of the previous section we can expand the opacity (k) to order e 2 /he, which results in a substantial 
simplification 
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where we have focussed on the propagation of the ip+ characteristics. The results for the ip- characteristics are 
identical in magnitude and follow from symmetry. 

From the form of Eq. we see that the opacity is zero for waves traveling with their magnetic field vectors 
perpendicular to the external field (_L mode). This result agrees with Bialynicka-Birula's analysis Q who found that 
although a wave in the _L mode readily generates waves in the || mode, a wave in the _L mode does not change to first 
order. These selection rules result from the CP invariance of QED and may be gleamed from the selection rules for 
photon splitting ||. 

Because our analysis tracks the evolution of only a single mode, we will calculate the opacity in the limit where the 
magnetic field of the wave is parallel to the external field. Waves in the || mode generate higher harmonics in the || 
mode but none in the _L mode. In this limit, 



B R = B and B F = 



and Eq. 42 simplifies further, 
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We define a dimensionless auxiliary function F(£) to characterize the opacity due to shocking 

n=-F(o~ = -miB 1 - 

c ts k 

We define Ib to be the characteristic length over which the magnetic field of the wave would change by B^. Iq is 
positive for sections of the wave where the magnetic field strength increases as its passes a stationary observer, or 
equivalently in the frame of the wave itself, Ib is negative in sections where the field strength decreases in the direction 
of propagation. 

The function F(£) may be expanded in the weak-field limit (£ < 0.5) yielding 
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where B n denotes the nth Bernoulli number. In the strong-field limit (£ > 0.5), we obtain 
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where C( x ) i s the Riemann Zeta function and we have used the expansions of Heyl and Hernquis t M . 

Note from Fig. || that F(£) is positive for all field strengths. Thus, from examination of Eq. we see that the 
opacity is positive (shocks will develop) in regions where the magnetic field is increasing toward the direction of 
propagation in the frame of the wave. F(£) also reaches a maximum near the critical field strength. 
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V. THE PHYSICAL SHOCK: JUMP CONDITIONS 



We expect a shock to develop when and where the value of the Riemann invariant is discontinuous. From Eq. [55] and 
using the result B = t±E, we see that the invariants are simply the electric and magnetic field strengths associated 
with the wave. Fig. [| schematically depicts the evolution of the wave. The shock begins to form at an optical depth 
of one where the field of the wave becomes discontinuous. 

As in fluid shocks, dissipative processes prevent the field strengths from becoming double valued. We use the 
Maxwell equal area prescription ( [Q) to calculate the shock profile after the characteristic analysis indicates that 
the field strengths become double valued. We start with a sinusodial wavefront, 

B(y, t) — — Bo sin(y — a + t) = —B sin vq (68) 

and obtain the following equation for the characteristics 

v(t) = vq + t sin vq (69) 

in the frame of the wave. Furthermore, for convenience we use the optical depth to shock formation as the time unit. 
From Fig. | and Eq. || we see that the wave evolves symmetrically about v — n. The position of the shock is given 



by the location which divides the double-valued regions into equal areas. By symmetry this occurs at v = tt. The 
wavefront at r = 2 is constructed in this manner. 

To determine the dissipation of energy by the shock, we calculate the mean power of the wave 

P=- cB 2 dv = -a+Bl / sin 2 v dv (70) 
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That is, the shock is located at v s — tt. For r < 1 the only real solution to Eq. 73 is vo, s = 7r. Therefore, before the 
shock forms the mean power in the wave is simply 1/2o+Bq. 

Unfortunately, in general this equation can only be solved numerically; however, two limits exist which can be 
treated analytically. As the shock just begins to develop vo ~ v near the shock, so vq, s ~ tt. If we expand Eq. |73| 
about uq,s = tt, we obtain 
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2 57T in 



(75) 



We find that as soon as the shock forms at r = 1, the shock begins to dissipate energy from the wave. Additionally, 
the dissipation does not begin abruptly. This first two terms in this expansion are accurate to ~ 1% for r — 1 < 0.2. 
At late times, we can find a solution to Eq. |7TJ such that v s « 0. Here we obtain 

v ,s ~ — ^-r (76) 

T + 1 

Ps 4ftif (77) 

This expression is accurate to one percent for r > 12. The upper panel of Fig. [5] depicts the energy dissipation soon 
after the shock forms. It is apparent that the dissipation begins smoothly. The lower panel shows the late evolution. 

Here, in the preceding analysis we have assumed that the field gradients are small both in our linearization and in 
our selection of the Heisenberg-Euler Lagrangian. Our linearization, specifically the assumption that the gradient of 
the fields is small relative to the fields breaks down when 
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(78) 



where A e is the Compton wavelength of the electron. The Heisenberg-Euler Lagrangian also breaks down when the 
field changes by Bk over scales similar to A e . In this limit, one must use more powerful techniques, such as Schwinger's 
proper-time method || to determine the effective Lagrangian. 

When the field changes dramatically over scales similar to A e , our fluid approximation breaks down. We estimate 
the thickness of the shock to be approximately A e . To further understand the properties of the shock, we derive the 
jump conditions across the shock discontinuity. 

We move to the rest-frame of the shock and insist on the continuity of the dual to the field tensor and the energy- 
momentum tensor: 



d v T^ = and <9„9^ = (79) 

where T pv = \e' iV ^F 5l {e.g. §) and 

The first condition follows from the gauge invariance of the fields. The second condition represents the conservation 
of energy and momentum. These jump conditions are equivalent to those used by Boillat Jlo| who insisted that the 
dual of field tensor be continuous and that the Euler-Lagrange condition be satisfied. 

For clarity, in contrast to the analysis of the preceding sections, we examine the energy-momentum and the field 
tensors of the combined wave and constant background field. Using the techniques outlined in |jTTf , we find the energy 
momentum tensor for non-linear electrodynamics. The canonical tensor (0 Miy ) is constructed from the Lagrangian by 
means of a Legendre transformation, 
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where A p is the potential four-vector of the electromagnetic field. We construct the more familiar symmetrized 
energy-momentum tensor (0 Miy ) by subtracting a total divergence, 
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where the final term of Eq. |83| is zero by the Euler-Lagrange condition. Using the definitions of / and K we obtain, 

6"" = A^F PP - 8J—F PP F v - g py C. (85) 
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where J = F^T pv . 

Simplifying this expression yields, 
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dl p \ dK J v ' 

which for C = — jF pa F pa , the linear case, yields 

qi* v = -^F pa F pa + F pp F p v . (87) 

in agreement with Itzykson and Zuber's result |yj ]. 

When determining shock jump conditions for a fluid, one calculates the velocity of the discontinuity relative to the 
rest frame of the fluid. In analogy to a fluid, we can associate the rest frame of the electromagnetic field with the 
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frame in which the energy-momentum tensor is diagonal. This frame exists if either of the two invariants (J and K) 
is non-zero. In this frame, the electric and magnetic fields are parallel and their magnitudes are given by 
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where the + sign is chosen for / > and — for / < 0. If we take the fields to point along the x-axis, we obtain 
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To apply this to the geometry of the previous sections (Fig. [I]), we will assume that K — and that the shock front 
is parallel to the x — z plane and traveling toward increasing y. 
The geometry leads to the jump conditions: 



[6 02 ] = 0, [9 22 ] = and [T 12 ] = 0. 



(92) 



across the shock. We calculate the components of the energy-momentum tensor in the shock frame through a boost. 
If K 7^ 0, the Lorentz transformation from the rest frames to the shock frame would include a rotation as well. By 
boosting the rest frames into the shock frame we get the following jump conditions 
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[e 22 ] = 
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v 2 ) 1/2 . The three 



where v is the speed that the diagonalizing frame is moving relative to the shock and 7 = (1 
conditions are physically conservation of energy and momentum flux, and the continuity of the electric field parallel 
to the surface of shock. Analogous jump conditions are given by p2| for relativistic fluid shocks. If we define, 



e = w — p — —£., p = C 
we obtain the following jump conditions 



2—1, w = -2—1 and n = -VI 
01 01 



[e 02 ] = [^vw] = 
[e 22 ] = [p + ~f 2 v 2 w] = 

[T 12 ] = [yvn] = 

Landau and Lifshitz p2] give the velocities of the diagonalizing frames relative to the discontinuity 
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and the equation of the shock adiabatic 
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where the subscripts 1 and 2 denote conditions on either side of the shock. The equation for the shock adiabatic is 
automatically satisfied to first order. 
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Taking Eq. 100 and Eq. 101 and assuming that the shock strength is a linear perturbation on the background field 
we get 
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where AI = I\ — Ii > 0, AS = £?i — £?2 > and we have used Eq. [36] to simplify the expression. 

It is more useful to find the speed at which the shocks travel relative to the external field. In general, the rest frame 
travels at a velocity that satisfies M 
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For K = 0, the rest frame of the electromagnetic field travels at 
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in the y-direction relative to external field behind and ahead of the shock, respectively. 

By summing the velocities of the rest frames relative to the external field and the shock relative to the rest frames, 
we find that the shock travels at the speed of light in the medium relative to the external field 



Wl,2 



(108) 



so we find in agreement with Boillat [|10| that the Heisenberg-Euler Lagrangian does not admit shocks to first order 
in the strength of the discontinuity, if shocks are strictly defined as discontinuities that travel at a speed other that 
the speed of "sound" in the medium. 



VI. CONCLUSIONS 



We have developed a relativistic fluid dynamic description of the electromagnetic field to derive the characteristic 
equations for electromagnetic waves in the presence of a strong external magnetic field. By using a analytic expressions 
for the effective Lagrangian of QED, we have obtained simple expressions to estimate the opacity of waves to shocking 
for arbitrary magnetic field strengths. Furthermore, by deriving a concise expression for the energy-momentum 
tensor for an arbitrary electromagnetic Lagrangian, we have calculated the shock jump conditions and find that the 
discontinuities travel at the speed of light in the medium relative to the external magnetic field. By extending this 
fluid dynamic description, we follow the eventual decay of a electromagnetic disturbance travelling through an intense 
magnetic field. 

For shocks to form from electromagnetic waves, not only is a strong external field required but also a source of 
coherent radiation. A prime location for electromagnetic shocks is the vicinity of a neutron star with field strengths 
approaching and exceeding the critical value and coherent electromagnetic Alfven waves. The study of the nonlinear 
corrections to the propagation of radiation through a plasma is beyond the scope of this work; however, we expect 
shock formation to be a hallmark of the nonlinear corrections of quantum electrodynamics and possibly an important 
process in the energy transmission near neutron stars. 
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Shocks form here 




FIG. 2. Characteristics and Shock Formation 



12 



B (G) 



1Q 13 1Q 14 10 15 10 16 




1 1 1 1 i i i i i 1 1 1 1 i i i i i 1 1 1 1 i i i i i 1 1 il i i i i i 1 1 1 

0.1 1 10 100 1000 

f 



FIG. 3. The figure depicts the auxiliary function as a function of £ 
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FIG. 4. Evolution of a wave toward shock formation in the 
comoving frame. The lower panel traces the power spectrum of 
at an optical depth of one-half, at an optical depth of one and at 
of!,- 2 . 



frame of the wave. The upper panel traces the wave in the 
the wave. The successive lines denote original wave, the wave 
an optical depth of two. The bold line shows a power spectrum 
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